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Abstract. The periodically varying Lorentz force of the periodic solar 
magnetic field generated by the solar dynamo can induce two kinds of mo- 
tions: torsional oscillations and periodic variations in the meridional circu- 
lation. Observational evidence now exists for both these kinds of motions. 
We discuss our ongoing efibrt in theoretically studying the variations of the 
meridional circulation. Then we present our theoretical model of torsional 
oscillations, which addresses the question why these oscillations start be- 
fore sunspot cycles at latitudes higher than where sunspots are seen. 
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1. Introduction 

The solar cycle is associated with a periodically varying magnetic field produced by 
the solar dynamo. This periodically varying magnetic field must give rise to a period- 
ically varying Lorentz force. We consider the possible kinds of motions induced by 
this periodically varying Lorentz force. 

To figure out the kinds of motions induced, we need to look at the Navier-Stokes 
equation with the periodically varying Lorentz force inserted in it. We can consider 
two kinds of motions separately: (i) the motions in the ^ direction, and (ii) the mo- 
tions in (r, 0) plane. It may be mentioned that the equation of continuity has to be 
satisfied along with the Navier-Stokes equation. The equation of continuity ensures 
that motions in r and 9 directions are coupled together. So we cannot treat those two 
directions separately. However, motions in the (p direction can be treated separately. 
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The basic motion in the cj) direction is the differential rotation, which can be mod- 
ified by the periodically varying Lorentz force. Such a periodic variation of the dif- 
ferential rotation — known as torsional oscillations — has been known for about three 
decades since its discovery at the surface by Howard & LaBonte (1980). On the other 
hand, the basic motion in the (r, 6) plane is the meridional circulation and only recently 
evidence has started coming that there is a periodic variation of the meridional circu- 
lation with the solar cycle (Hathaway & Rightmire 2010). Soon after the discovery 
of torsional oscillations, some authors started exploring their theoretical impUcations 
(Yoshimura 1981; Schiissler 1981; Tuominen & Virtanen 1984). Due to the rapid 
developments in dynamo theory in the last few years, this subject needs revisiting. 

We (Karak & Choudhuri) are working on a theoretical model of the variation 
of meridional circulation due to the Lorentz force. Some of the basic concerns are 
discussed in §2. Then in §3 we suimnarize our recent work on torsional oscillations 
(Chakraborty, Choudhuri & Chatterjee 2009a). 

2. Variations of meridional circulation 

The meridional circulation of the Sun plays a crucial role in the flux transport dynamo. 
However, its theory is still rather poorly understood. HeUoseismology provides some 
information about the meridional circulation in the upper layers of the convection 
zone (Giles et al. 1997; Braun & Fan 1998; Gonzalez Hermandez et al. 2006; Svanda, 
Kosovichev & Zhao 2007). But there is no reUable observational data available yet on 
the nature of the meridional circulation at the bottom of the convection zone, which 
is crucial for the dynamo (Nandy & Choudhuri 2002). 

The classic investigation by Kitchatinov & Rudiger (1995) showed that the meri- 
dional circulation arises out of a slight imbalance between two large terms, requiring 
a pole-equator temperature diflFerence of about 5°. A consequence of this is that even 
small fluctuations in any of these large terms may cause significant variations in me- 
ridional circulation. On the basis of indirect evidence, Karak (2010) concluded that 
there have been large variations in meridional circulation in the last few centuries. 
The variations in meridional circulation also seem to be the main reason behind the 
well-known Waldmeier effect of solar cycles (Karak & Choudhuri 201 la). 

Apart from these random variations of meridional circulation, we expect some 
systematic periodic variations with the solar cycle. Only recently Hathaway & Right- 
mire (2010) and Basu & Antia (2010) have presented evidence for this. It appears 
that the meridional circulation near the surface becomes somewhat weaker at the time 
of the sunspot maximum. We expect a strong toroidal field at the bottom of the con- 
vection zone at the time of the sunspot maximum. Such a toroidal field has a Lorentz 
force directed towards the rotation axis and has a tendency of slipping poleward (van 
Ballegooijen & Choudhuri 1988). This poleward shpping tendency will oppose the 
equatorward meridional circulation at the bottom of the convection zone. We (Karak 
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& Choudhuri) are now carrying out a detailed calculation to investigate whether the 
surface observations of the reduction of meridional circulation at the time of sunspot 
maximum can be explained as arising out of this opposition to meridional circulation 
due to the poleward slip tendency of the strong toroidal field at the bottom of the 
convection zone. 

It may be mentioned that Nandy, Munoz-Jaramillo & Martens (2011) have as- 
sumed in a recent work that the meridional circulation changes randomly at every 
solar maximum. We disagree with this assumption and believe that the meridional 
circulation decreases at the solar maximum due to the Lorentz force of the magnetic 
fields in a systematic deterministic way. If the Lorentz force quenches the meridional 
circulation at the time of the solar maximum, one important question is whether this 
will have any efi'ect on the dynamo. Recently Karak & Choudhuri (201 lb) performed 
a dynamo simulation by assuming that the ampUtude of the meridional circulation is 
quenched according to the equation 

i;o = v'JW + (B/Bofl (I) 

where B is the average value of the toroidal field at the bottom of the convection zone. 
If the value of the turbulent diffusivity is assumed to be sufficiently small (rj ~ 10'"- 
lO" cm^ s '), then the dynamo becomes unstable on including such a quenching of 
the meridional circulation. The dynamo remains stable only if the turbulent diffusivity 
is in the range ~ lO'^-lO'^ cm^ s"' . This provides another argument in favour of a 
high-diffusivity dynamo, strengthening a case already made by several authors (Jiang, 
Chatterjee & Choudhuri 2007; Goel & Choudhuri 2009; Hotta & Yokoyama 2010; 
Karak 2010). 

3. Torsional oscillations 

The small periodic variation in the Sun's rotation with the sunspot cycle, first dis- 
covered on the solar surface by Howard & LaBonte (1980), is called torsional oscil- 
lations. Helioseismology has now established its existence throughout the convection 
zone (see Howe et al. 2005 and references therein). Its amplitude near the surface is of 
order 5ms' or about 1 % of the angular velocity. Apart from the equatorward-propa- 
gating branch which moves with the sunspot belt, there is also a poleward-propagating 
branch at high latitudes. One intriguing aspect of the equatorward-propagating branch 
is that it begins a couple of years before the sunspots of a particular cycle appear and at 
a latitude higher than where the first sunspots are seen. The top panel of Fig. 1 shows 
the torsional oscillations at the solar surface with the butterfly diagram of sunspots. 
If the torsional oscillation is caused by the Lorentz force of the dynamo-generated 
magnetic field as generally believed, then the early initiation of this oscillation at a 
higher latitude does look like a violation of causality! The main aim of our recent 
work (Chakraborty, Choudhuri & Chatterjee 2009a, 2009b) has been to explain this 
which could not be explained by the earlier theoretical models (Dumey 2000; Covas 
et al. 2000; Bushby 2006; Rempel 2006). 
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Figure 1. The time-latitude plot of torsional oscillations on the solar surface with the buttertly 
diagram of sunspots superposed on it. The upper panel is based on observational data of surface 
velocity measured at Mount Wilson Observatory (courtesy: Roger Ulrich). The bottom panel 
is from our theoretical simulation. 

Our calculations are based on the dynamo model presented by Chatterjee, Nandy 
& Choudhuri (2004). In order to model torsional oscillations, in addition to the basic 
equations of the dynamo, we simultaneously have to solve the <p component of the 
Navier-Stokes equation in the form 



P&+£).M^ = £)vK] + (Fz.) 



(2) 



where D„[v^] is the term corresponding to advection by the meridional circulation, 
is the diffusion term, and (Fi)^ is the component of the Lorentz force. If 
the magnetic field is assumed to have the standard form 



B = Bir, 0, f)e^ + V x [A(r, 0, f)e^ 
then the Lorentz force is given by the Jacobian 



47r(Fi)^ = -7 



1 lsB4„sA 



r,0 



(3) 



(4) 



where s = rsinfl. On the basis of flux tube simulations suggesting that the mag- 
netic field in the tachocline should be of order 10^ G (Choudhuri & Oilman 1987; 
Choudhuri 1989; D'Silva & Choudhuri 1993), it is argued by Choudhuri (2003) that 
the magnetic field has to be intermittent in the tachocline. Hence the full expression 
of the Lorentz force involves a filling factor as explained by Chakraborty, Choudhuri 
& Chatterjee (2009a). 
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Figure 2. The depth-time plot of torsional oscillations at latitude 20°. The left panel 
from Vorontsov et al. (2002) is based on SOHO observations, whereas the right panel from 
Chakraborty, Choudhuri & Chatterjee (2009b) is based on our theoretical simulation. The solid 
and dashed lines in the right panel indicate the Lorentz force (positive and negative values 
respectively). 



Our theoretical model incorporates a hypothesis proposed by Nandy & Choudhuri 
(2002), which is essential for explaining the early initiation of the torsional oscillation 
at high latitudes. According to this Nandy-Choudhuri (NC) hypothesis, the meridi- 
onal flow penetrates into the stable layers below the convection zone at high latitudes. 
This causes the formation of toroidal field in the high-latitude tachocline. Sunspots 
form a few years later when this field is advected to lower latitudes and brought in- 
side the convection zone. We also assume that the stress of the magnetic field formed 
in the tachocline is carried upward by Alfven waves propagating along vertical flux 
concentrations conjectured by Choudhuri (2003). 

The incorporation of the NC hypothesis in our theoretical model causes magnetic 
stresses to build up at higher latitudes before sunspots of the cycle appear, leading to 
the early initiation of torsional oscillations. The bottom panel of Fig. 1 shows the- 
oretical results of torsional oscillations at the surface with the theoretical butterfly 
diagram. This bottom panel can be compared with the observational upper panel in 
Fig. 1. Our theoretical model also gives a satisfactory account of the evolution of 
torsional oscillations within the convection zone. The depth-time plot of torsional 
oscillations at a certain latitude given in Fig. 3 of Chakraborty, Choudhuri & Chat- 
terjee (2009b) compares favourably with the observational plot given in Fig. 3(D) of 
Vorontsov et al. (2002). This is reproduced here in Fig. 2 for completeness. We end 
by pointing out that our theoretical model of torsional oscillations attempts to ex- 
plain various aspects of observational data in much greater detail than what had been 
attempted in previous theoretical models. 
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